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The motion of an optically trapped sphere constrained by the vicinity of a wall is investigated 
at times where hydrodynamic memory is significant. First, we quantify, in bulk, the infiuence of 
confinement arising from the trapping potential on the sphere's velocity autocorrelation function 
C{t). Next, we study the splitting of C(t) into C||(t) and C±{t), when the sphere is approached 
towards a surface. Thereby, we monitor the crossover from a slow t~^^^ long-time tail, away from 
the wall, to a faster t~^^'^ decay, due to the subtle interplay between hydrodynamic backfiow and 
wall effects. Finally, we discuss the resulting asymmetric time-dependent diffusion coefficients. 
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Understanding and controlling the transport of col- 
loidal microcarriers, such as membrane vesicles, through 
a fluid is one of the main challenges in cell biology [H , and 
related lab-on-a-chip approaches. With the miniaturiza- 
tion of such technologies, particularly microfluidics 
the colloid's motion is increasingly confined, and the in- 
fluence of boundaries, e.g. a channel wall, becomes non- 
negligible. In particular, in a system as small as a cell, 
many obstacles will alter the trajectory of a diffusing par- 
ticle. Any deviation from its well-understood free Brow- 
nian motion, will give information on the particle's sur- 
roundings. The reduction of a colloid's mobility close to 
a wall, also known as 'surface conflnement', was already 
predicted by Lorentz in 1907 0, 13|, and is expected to 
entail drastic effects on its time-dependent Brownian mo- 
tion arising from the thermal fluctuations in the system. 

Experimental evidences for this wall effect are compa- 
rably recent, and consist mostly in measuring changes 
in the diffusion coefficient I? of a micron-sized sphere 
approaching a surface [5, ,6,, or being confined between 
two walls [81]. Measurements confirm that an interface in- 
creases the steady-state friction, hence the viscous drag 
on the particle, slowing down its diffusion as the dis- 
tance h between the sphere's center and the surface is 
reduced. The motion becomes anisotropic in the direc- 
tions parallel and perpendicular to the wall. Conven- 
tional data acquisition and analysis rely on Lorentz 's 
zero-frequency approximation for the mobility, which 
only takes a viscous contribution from the fluid into ac- 
count. However, for a neutrally buoyant sphere, the 
presence of hydrodynamic memory due to momentum 
conservation in the fluid [1, El leads to it having 
a signiflcantly delayed dynamic behavior, at time scales 
much larger than the particle's momentum relaxation 
time Tp = 2a^pp/977. This memory has been observed 
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in bulk for colloidal suspensions [12 



cently, directly for a single microsphere 




lies in the back-flow that a spherical particle of density 
Pp and radius a creates in a fluid of viscosity rj and den- 
sity P(. The long time diffusion of momentum in the 
viscous fluid leads to an algebraic decay, the so-called 
long-time tail, in the velocity autocorrelation function 
(VACF) Cit) = {v{t)v{0)) ~ L»y^i74^t-3/2 of the fluc- 
tuating sphere at times t >T{ ~ pfO^/r], the time needed 
by the perturbed fluid flow fleld to diffuse over the dis- 
tance of one particle radius. Obviously, the presence of 
a wall bounding the fluid will affect the fluid vortex once 
it has encountered the wall, which will occur at times 
larger than its propagation time, r„ = h^pf/rj, from the 
particle to the boundary [l^ . 

In this Letter, we investigate the effect of hydrody- 
namic vortex diffusion on the nature of the long-time tail 
as a function of the wall-particle distance h, by measur- 
ing the position fluctuations of a single spherical particle 
approaching a hard surface. Weak optical trapping [l3| 
is employed to position a silica sphere relative to a plain 
glass surface, and, at the same time, to track the Brow- 
nian particle's trajectory interferometrically [l3| with a 
precision better than 1 nm and a time resolution of 2 jis 
[1J|. The optical trap is created by focusing a 20 x- 
expanded Nd-YAG beam (A = 1064 nm) by a 63 x water- 
immersion objective lens (NA= 1.2). The thermal posi- 
tion fluctuations of the sphere in the weak trap are mea- 
sured with an InGaAs quadrant photodiode, amplifled 
and digitized (12 bits). The position signal is recorded 
during is = 20 s with a sampling rate fs = 500 kHz corre- 
sponding to iV = 10'' data points. The Brownian particle 
used here is a silica sphere of several micrometers diam- 
eter {pp — 1.96 g/cm"^) immersed in water (pf = Ig/cm"^, 
?7 = 10^'^Pa-s). It is gradually approached towards the 
surface of a 100 jim-sized sphere sandwiched between the 
two coverslides of a fluid chamber (size « 2 cm x 0.5 cm 
and thickness « 100 pm). As illustrated in Fig.[Tl the big 
sphere can be considered as a flat surface on the scale of 
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FIG. 1: (color) (a) 3D lateral view of the experiment; spheres 
are drawn to scale. A silica particle of radius a =1.5 pm 
trapped by the laser focus is placed next to the surface of 
a significantly larger silica sphere. This 100 pm sphere is im- 
mobilized between the two coverglass surfaces of the sample 
chamber, (b) Optical image of the probing particle's position 
relative to the wall created by the big sphere. The 3 pm prob- 
ing particle was placed at a distance h =11.5 pm away from 
the 100 pm sphere's surface and gradually approached. The 
velocity correlation functions, as well as the diffusion coeffi- 
cients for the motion parallel and perpendicular to the wall 
are measured. 



our Brownian sphere. Such a configuration circumvents 
the draw-back of the lower resolution along the optical 
axis z, intrinsic to optical trapping interferometry [l9| . 
thereby allowing measurements with comparable accu- 
racy in both directions, parallel and perpendicular to 
the boundary. The small silica spheres are suspended 
in water at a low enough concentration to allow trapping 
and observation of one isolated particle. The sample is 
mounted onto a piezo-stage, and the 100 pm sphere can 
be positioned at a distance h relative to the trapped par- 
ticle by moving the piezo-stage in all three dimensions 
with a precision of w 1 nm. 

As a first test of the sensitivity of our set-up, we 
measure the VACF for a microsphere (a = 2.25 jim, 
Tp ~ 2.2 ps, Tf = 5.1 ps) far away from any boundaries 
{h >40pm). The optical trap is particularly well-suited 
for the experiment, as its harmonic potential results in 
a linear force F ~ —kx{t), which only acts on the par- 
ticle, but does not confine the surrounding fiuid. The 
spring constant k of the trap gives rise to a new time scale 
Tk = Q-K-qa/k. A model for the Brownian motion of a par- 
ticle in a harmonic potential is provided by Clercx and 
Schramm [iOl , who went beyond a simple Stokes approx- 
imation by using the time-dependant linearized Navier- 
Stokes equation to describe the fiuid motion. Fig.Oshows 
the measured normalized velocity autocorrelation func- 
tion C(i)/C(0) for k = 19pN/m and k = 5pN/m. In 
order to compare data with the theoretical curves, the 
normalization C(0) and k are treated as fit parameters. 
As can be seen in Fig. [21 our data follow the theoretical 
prediction [i^l over the full time range, down to a noise 
level of 5 • 10~^. As was shown earlier jS^I, an inter- 
mediate time window Tf < t < T]^/20 opens, where the 




FIG. 2: (color) Log-log plot for the measured normalized ve- 
locity autocorrelation of a sphere (a = 2.25pm) held by an 
optical trapping potential with k = 19pN/m (green circles) 
and with k = 5pN/m (red squares). The continuous red and 
green lines are fits to the model of Clercx and Schramm (20l ]. 
The hydrodynamic long-time tail C{t) oc t~^^^ (black line) 
emerges as the spring constant is decreased. Inset: Blow-up 
of the region of anti-correlations on a log-linear scale. 



particle's motion can still be observed, but is free from 
the infiuence of the trap and mainly dominated by the 
fiuid's inertia. Hence, in order to minimize infiuences 
from the trapping potential on the Brownian motion and 
resolve hydrodynamic memory effects, we adjust Tk to be 
as long as possible. By decreasing the spring constant 
from k = 19pN/m to k ~ 5pN/m, we directly observe 
the emergence of the long-time anomaly t~^^'^ (Fig- E]) 
over two decades in signal, which has not been previ- 
ously reported in the literature. In an incompressible 
liquid, the initial value of the VACF is determined by 
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where the amplitude B = (pp/pf + l/2)/9^/7r depends 
only on the mass ratio. The value of the prefactor is of 
the order B ~ 1/10 for the colloidal system used here, 
and can be directly read off the data at t = Tf (Fig. [2]) . 
At longer times, t 3> Tk (t = 22ms for k = 19pN/m 
and t = 96ms for k = 5pN/m), the fiuid's inertia is ex- 
pected to generate a second zero in C (t) followed by a fast 
t-^/2-tail as (15S/4)(t/Tf)-3/2(t/rk)-2^ which, however 
remains unobservable due to noise. 

Next, we approach the Brownian sphere (a = 1.5 pm, 
Tp = Ips, Tf = 2.25 ps) towards the boundary created 
by the 100pm glass sphere (Fig. [J), and vary h from 
37.8 pm to 4.8 pm, corresponding to a reduction in t„ 
from 1400 ps to 23 ps. The trap stiffness is minimized 
to fc « 2pN/m (Tk =14ms) by lowering the incoming 
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laser power. It is worth noting that decreasing k de- 
grades the signal-to-noise level, as can be seen in the 
increase of the error bars in Fig. [51 As already men- 
tioned above, momentum from the fluid is transferred to 
the wall at times t > t^. The leading hydrodynamic 
tail t"'^/^ in C{t) is then canceled, and C{t) splits into 
C|| (t) for the motion parallel and C± {t) for the motion 
perpendicular to the wall. Recently, Felderhof pro- 
vided a full analytical solution for the motion in both 
directions. He generalized the frequency-dependent ad- 
mittance {i.e. the frequency-dependent mobility) for the 
unconstrained motion 3^o('^) [Ul to 3^|[(w) and y±{uj) 
for a point particle moving close to a wall. Relying on 
the fluctuation-dissipation theorem j2l|, he then calcu- 
lated C||(t) and C±{t) via a Fourier-back-transform from 
3^11 (w) and y±{uj). For the parallel motion, a more rapid 
but still algebraic decay is predicted 
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for t > provided that the wall is not too close; r„ » 
Tf, Tp. For the motion perpendicular to the wall, the long- 
time behavior, t > r^, is predicted, to leading order in 

a/h, as 
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where A^s = {2p^/pi - 5)/9. 

To disentangle surface confinement from the trap con- 
straint, we extend Felderhof's theoretical app roach 
by including the optical trapping force 12311 . In bulk, 
the harmonic potential modifies 3^o(w) to 3^k(w) = 
[3^o(^)~^ — k/\ijj]~^ jiOl- Accordingly, in the presence of 
a wall, the admittances 3^||,_L(a;) have to be modified to 
[3^11 _L(a')~"'^ — k/\ijj\~^, and a suitable Fourier algorithm 
yields the VACFs in the entire time domain, which is 
fitted to our data. 

Figure [3] shows the normalized C||(i) and C^(t), in the 
time range Tf < i < ti./20 = 0.7 ms, at four different 
distances from the wall. Far away from the boundary, 
at h =37.8 pm, we observe the same t~'^/^ power-law for 
both directions (Fig.[3K). However, as soon as the sphere 
reaches the proximity of the wall, and r„ falls into our 
window of observation, C{t) splits due to the hydrody- 
namic interaction with the wall (Fig. [3b, h = 9.8 ym, and 
T„ = 96 ys). Our data show clearly that the anisotropy 
increases even more as h is further decreased (Fig. [3fc 
and d). The higher noise floor at smaller h probably 
results from scattering of the highly divergent trapping 
laser beam by the big sphere. In C||(t), a transition 
from the free bulk behavior at h =37.8 pm, character- 
ized by the t"'^/^ power-law (Fig. [SK), to confined mo- 
tion, with a steeper, power-law arises in the data 
at t ~ Tto. The effective amplitude in Eq. ([2]) is reduced 
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FIG. 3: Log-log plot of both normalized VACF, (t)/C|| (0) 
and C±{t)/C±{0) for a sphere (a =1.5 pm, Tp = 1 ps, T{ = 
2.25 ps) trapped in a weak optical potential {k ^ 2pN/m, 
Tk = 14 ms ). The increasingly anisotropic VACF is measured 
at four distances from the wall (/i=37.8, 9.8, 6.8 and 4.8 pm, 
corresponding to r„ = 1400, 96, 46, 23 ps, respectively). The 
characteristic power-laws from Eq. ^ and Eq. ([S} are repre- 
sented by black lines as guides to the eye. The experimental 
data (squares and circles) are compared to the theory that 
includes hydrodynamic memory, wall effects, and harmonic 
restoring forces (continuous lines). 



by up to 30% due to the trap. The drop to negative val- 
ues in C||(i) arises from anticorrelations imposed by the 
harmonic trapping potential, and can only be captured 
by our extended theory As can be inferred from 

Figs. [HI the theoretical curves describe the data down to 
the noise level. For the perpendicular motion, the pref- 
actor Az3 of the leading f~^/^-term in Eq. ([3]) depends 
on the relative densities, and is negative in our experi- 
ment, Az3 = —0.12. Hence, for k = 0, a sign change 
is expected at t ~ -tI/4t{A^3 = 1890,8.53,1.96, and 
0.49 ms. According to Eq. ([3]), a positive tail i"''/^ should 
dominate in the time window t„ < t <C t'^/t{. How- 
ever, the sign change in C\^{t) is observed to occur much 
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FIG. 4: Normalized time-dependent diffusion coefficients, 
D\\{t)/D and D±{t)/D, for the direction parallel and perpen- 
dicular to the wall at the same distances h as in Fig. [S] The 
arrows at the right correspond to the asymptotic values given 
by Lorentz's prediction in Eq. ([4]). The experimental data are 
represented by symbols, whereas the full lines correspond to 
the theoretical fits. The arrows at the top and the bottom 
indicate the respective time-points of the shallow maxima. 



earlier, even though the trapping potential is minimal. 
Nevertheless, for the two closest distances (Fig. [S]: and 
d) , the steeper decay is compatible with the intermediate 
power law i"^/^ predicted by Eq. ([3]), but with an am- 
plitude reduced to 50% by the trap. This influence from 
the trap as well as noise also obscure the crossover from 
the positive to the negative tail. 

In general the flts in Fig. [3] are not very sensitive 
to the spring constant k. Therefore, we determine its 
value by comparing our data to the long-time behavior 
of the time-dependent diffusion coefficients, D^\^j_{t) = 

J^Cij_{t')dt', shown in Fig. H As can be inferred from 
Fig. m the suppression of diffusion becomes anisotropic 
and follows the point-particle prediction by Lorentz [3|, 

Z?!! = D[l - 9a/16h] , Di_ = D[l - 9a/Sh] , (4) 

where D = ksT /Qirrja is the diffusion constant in bulk. 
The higher orders in a/h are known to contribute less 
than 2% for /i/a > 3 y, and our data analysis suggests 
that the point-particle Hmit is also valid for the time- 
dependent motion with similar accuracy. The observed 
zero-crossings in the VACF translate into shallow max- 
ima in the diffusion coefficients. They are still up to 5% 
below the asymptotic values D^^ and D±, expected for 
/c = 0, which exemplifles that free Brownian motion is 
not attained in optical trapping of micron-sized spheres, 
even for times t < T]^. At longer times the trap domi- 
nates and reduces the diffusion coefHcients to zero (not 
shown). 

In conclusion, we provide high precision experimen- 
tal data relying on optical trapping interferometry which 
validate recent theoretical models. This approach gives 



access to the study of diffusion effects at a boundary, 
which are particularly relevant in diffusion-mediated in- 
teractions. During such processes, structures, like a large 
protein or membrane, interact with each other. This 
search for a potential interaction partner may be favored 
by the wall effect described here, as a fairly high mobility 
is preserved along the boundary, whereas diffusion away 
from the interaction partner is suppressed [23 |. 
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